The role of exponential asymptotics and complex singularities in self-similarity, transitions, and branch merging of nonlinear dynamics

We study a prototypical example in nonlinear dynamics where transition to self-similarity in a singular limit is fundamentally changed as a parameter is varied. Here, we focus on the complicated dynamics that occur in a generalised unstable thin-film equation that yields finite-time rupture. A parameter, n, is introduced to model more general disjoining pressures. For the standard case of van der Waals intermolecular forces, n = 3, it was previously established that a countably infinite number of self-similar solutions exist leading to rupture. Each solution can be indexed by a parameter, ϵ = ϵ1 > ϵ2 > · · · > 0, and the prediction of the discrete set of solutions requires examination of terms beyond-all-orders in ϵ. However, recent numerical results have demonstrated the surprising complexity that exists for general values of n. In particular, the bifurcation structure of self-similar solutions now exhibits branch merging as n is varied. In this work, we shall present key ideas of how branch merging can be interpreted via exponential asymptotics

Anomalous diffusion in polymers: long-time behaviour

We study the Dirichlet boundary value problem for viscoelastic diffusion in polymers. We show that its weak solutions generate a dissipative semiflow. We construct the minimal trajectory attractor and the global attractor for this problem.Comment: 13 page

Dynamics and universal scaling law in geometrically-controlled sessile drop evaporation

The evaporation of a liquid drop on a solid substrate is a remarkably common phenomenon. Yet, the complexity of the underlying mechanisms has constrained previous studies to sphericallysymmetric configurations. Here we investigate well-defined, non-spherical evaporating drops of pure liquids and binary mixtures. We deduce a universal scaling law for the evaporation rate valid for any shape and demonstrate that more curved regions lead to preferential localized depositions in particle-laden drops. Furthermore, geometry induces well-defined flow structures within the drop that change according to the driving mechanism. In the case of binary mixtures, geometry dictates the spatial segregation of the more volatile component as it is depleted. Our results suggest that the drop geometry can be exploited to prescribe the particle deposition and evaporative dynamics of pure drops and the mixing characteristics of multicomponent drops, which may be of interest to a wide range of industrial and scientific applications

The role of Allee effect in modelling post resection recurrence of glioblastoma

Resection of the bulk of a tumour often cannot eliminate all cancer cells, due to their infiltration into the surrounding healthy tissue. This may lead to recurrence of the tumour at a later time. We use a reaction-diffusion equation based model of tumour growth to investigate how the invasion front is delayed by resection, and how this depends on the density and behaviour of the remaining cancer cells. We show that the delay time is highly sensitive to qualitative details of the proliferation dynamics of the cancer cell population. The typically assumed logistic type proliferation leads to unrealistic results, predicting immediate recurrence. We find that in glioblastoma cell cultures the cell proliferation rate is an increasing function of the density at small cell densities. Our analysis suggests that cooperative behaviour of cancer cells, analogous to the Allee effect in ecology, can play a critical role in determining the time until tumour recurrence

Equilibrium interface solutions of a degenerate singular Cahn-Hilliard equation

AbstractWe present an analysis of the equilibrium diffusive interfaces in a model for the interaction of layers of pure polymers. The discussion focuses on the important qualitative features of the solutions of the nonlinear singular Cahn-Hilliard equation with degenerate mobility for the Flory-Huggins-de Gennes free energy model. The spatial structure of possible equilibrium phase separated solutions are found. Using phase plane analysis, we obtain heteroclinic and homoclinic degenerate weak compact-support solutions that are relevant to finite domain boundary value problems and localized impurities in infinite layers

Intermediate asymptotics for Richards' equation in a finite layer

Perturbation methods are applied to study an initial-boundary-value problem for Richards' equation, describing vertical infiltration of water into a finite layer of soil. This problem for the degenerate diffusion equation with convection and Dirichlet/Robin boundary conditions exhibits several different regimes of behavior. Boundary-layer analysis and short-time asymptotics are used to describe the structure of similarity solutions, traveling waves, and other solution states and the transitions connecting these different intermediate asymptotic regimes

Similarity solutions of the lubrication equation

AbstractWe present a method for constructing closed-form similarity solutions of the fourth-order nonlinear lubrication equation. By extending a technique used to study second-order degenerate diffusion problems, corresponding interface profiles and diffusion coefficient functions can be derived in exact form. Different classes of spreading and shrinking solutions are obtained using this approach

Motion of wetting fronts moving into partially pre-wet soil

We study the motion of wetting fronts for vertical infiltration problems as modeled by Richards' equation. Parlange and others have shown that wetting fronts in infiltration flows can be described by traveling wave solutions. If the soil layer is not initially dry, but has an initial distribution of water content then the motion of the wetting front will change due to the interaction of the infiltrating flow with the pre-existing soil conditions. Using traveling wave profiles, we construct simple approximate solutions of initial-boundary value problems for Richards' equation that accurately describe the position and moisture distribution of the wetting front. We show that the influences of surface boundary conditions and initial conditions produce shifts to the position of the wetting front. The shifts can be calculated by examining the cumulative infiltration, and are validated numerically for several problems for Richards' equation and the linear advection-diffusion equation. © 2005 Elsevier Ltd. All rights reserved